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ABSTRACT 

In this paper we use recent X-ray and radio observations of the ejecta-rich Galactic supernova 
remnant (SNR) G337.2— 0.7 to determine properties of the supernova (SN) explosion that formed this 
source. H I absorption measurements from the Australia Telescope Compact Array (ATCA) constrain 
the distance to G337.2— 0.7 to lie between 2.0 ± 0.5 and 9.3 ± 0.3 kpc. Combined with a clear radio 
image of the outer blast-wave, this distance allows us to estimate the dynamical age (between 750 
and 3500 years) from the global X-ray spectrum obtained with the XMM-Newton and Chandra 
observatories. The presence of ejecta is confirmed by the pattern of fitted relative abundances, which 
show Mg, Ar and Fe to be less enriched (compared to solar) than Si, S or Ca, and the ratio of Ca 
to Si to be 3.4±0.8 times the solar value (under the assumption of a single electron temperature and 
single ionization timescale). With the addition of a solar abundance component for emission from the 
blast-wave, these abundances (with the exception of Fe) resemble the ejecta of a Type la, rather than 
core-collapse, SN. Comparing directly to models of the ejecta and blast- wave X-ray emission calculated 
by evolving realistic SN la explosions to the remnant stage allows us to deduce that one-dimensional 
delayed detonation and pulsed delayed detonation models can indeed reproduce the major features of 
the global spectrum. In particular, stratification of the ejecta, with the Fe shocked most recently, is 
required to explain the lack of prominent Fe-K emission. 

Subject headings: supernova remnants: individual (SNR G337.2— 0.7) — supernovae: general — nu- 
clear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

A supernova (SN) explosion marks the death of a star, 
whether that be the gravitational collapse of a massive 
star that has run out of nuclear fuel at its core (core- 
collapse) or the disruption of a white dwarf by runaway 
nuclear reactions (Type la). The long-lived remnants 
of SN explosions provide a rich tapestry of information 
about the metal-rich ejecta produced in a SN, as well 
as the environment into which the remnant is expand- 
ing. X-ray spectra of supernova remnants (SNRs) are 
particularly useful given the extremely high shock tem- 
peratures and the fact that most cosmically abundant 
elements produced in the explosion have diagnostic K- 
shell lines in the 0.5 to 8.0 keV X-ray bandpass. Fur- 
thermore, the expansion of the ejecta allows their struc- 
ture to be spatially resolved. Because X-ray emission 
from the ejecta can be identified thousands of years af- 
ter the explosion, there exists a sample of young nearby 
SNRs despite a galactic SN explosion rate of only ~1 per 
century Given the diversity within the classes of core- 
collapse and Type la explosions, and the variety of open 
issues regarding their explosion mechanisms, it behooves 
us to study all the young SNRs in our Galaxy and nearby 
in the hope of gaining an understanding of the supernova 
explosions that originated them. 

Evidence for metal-rich ejecta in SNR G337.2— 0.7 was 
first found based on the strength of the emission lines 
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in the ASCA X-ray spectrum (Rakowski ct al. 2001). 
In that analysis, Si and S were confidently shown to be 
over-abundant compared to their solar abundances, and 
the pattern of metal abundances relative to each other 
was also highly non-solar. However the abundance ratios 
between the metals were not sufficiently well constrained 
to determine the origin of the SNR. In this paper, we 
analyze new observations of G337.2— 0.7 with Chandra 
and XMM-Newton, and we complement these X-ray ob- 
servations with radio data from the Australia Telescope 
Compact Array (ATCA). In our analysis, we devote spe- 
cial attention to the limitations of the techniques that are 
commonly used to measure elemental abundances, and 
how they influence our ultimate goal of determining the 
origin of this supernova remnant and contributing to our 
current understanding of SN explosions. Great care must 
be taken when interpreting X-ray spectra from SNRs be- 
cause this emission will be dominated by the hottest, 
densest highly ionized material. Abundant elements can 
be hidden if they are in a low density low temperature 
region, not co-mingled with the bulk of the ejecta. Even 
amongst the bright ejecta, if the standard single electron 
temperature, single ionization timescale models are used, 
hot plasma may be misinterpreted as overabundances. 

From the theoretical side, numerical simulations of SNe 
are rapidly increasing in sophistication. The inclusion of 
magnetic fields, angular momentum, asymmetries, and 
instabilities in the flow all require a 3-D treatment and 
are beginning to be explored in both core-collapse and 
Type la explosions. In this introduction we highlight 
some of the core issues under debate in the theory of SN 
explosions, which will set the context for our in-depth 
analysis of G337.2-0.7. 



A key question for core-collapse explosions is what re- 
juvenates the shock after core-bounce. When the in- 
falling material forces the core of the star to such a 
high density that a neutron-rich nuclear composition 
is favored, the core stiffens and pushes back ("core- 
bounce"). However, in spherically symmetric simula- 
tions, this shock will stall against the force of the con- 
tinually infalling material and fail to cause an explo- 
sion. A multi-dimensional mechanism is required to 
revitalize the shoc k, with the co ntenders includi ng jet- 
induced explosions (Khokhlqv_et^Ljy.999; W heeler et al.l 
l207r3lMa,eda, et al.ll2002HAkiva,ma et a,l]l2003ft. neutrino 
driven convection ilFrver fc Hegerl 120001 L Kifonidi s et al.l 
l2003t iKotake et all 12003: Bur rows et alJ 12005ft . and 
turbulent dissipation of magneto-rotational energy 
^Thompson et all [2005). Magneto- hydrodynamically 
(MHD) driven jets have particularly been proposed as 
a mechanism for gamma-ray burst (G RB) explosions 
(jWheeler et al.ll2002HMaeda et al.ll2002ft . The magnetic 
field may be pre-existing or it c ould be amplified b y 
the magneto-rotational instability l)Akivama et alJ2 003). 
There have also been improvement s in neutrino trans- 
port calculations IjLivne et al.ll2004ft and the incorpora- 
tion of rotation into neutrino-driven convection models 
IlFrver fe Hegerll2000l IKotake et alJl2003l iBurrows et all 
2005). Including rotation naturally leads to a bi-polar 
structure in the neutrino flux and accretion of mass, such 
that the increased neutrino heating is sufficient and di- 
rected enough to revitalize the shock near the poles and 
explode the star. Alternatively, the rotational energy 
of the in-falling material itself can be converted into an- 
other source of heating behind the stalled shock by turbu- 
lent dissipation of the differential shear via th e magneto- 
rotational instability ({Thompso n et al.ll2005ft . 

Simulations of white dwarf explosions are also pro- 
gressing. The main difficulty is in constructing a realistic 
explosion that reproduces the abundances and stratifi- 
cation of the ejecta inferred from Type la spectra (e.g. 
iBranch et al.l 12005. and references therein). In one di- 
mension, supersonic flame propagation (detonation) can 
burn almost the entire star to Fe-group elements be- 
cause the star retains the original density ahead of the 
flame. A purely sub-sonic flame (deflagration) will push 
the star outwards, reducing the density of the mate- 
rial ahead of the shock. This causes the flame to be 
quenched earlier than in the detonation case and leave 
behind a large outer shell of unburnt material. 4 Neither 
pure detonation nor pure deflagration produces the nec- 
essary quantities of intermediate mass elements (between 
O and Fe). The ad hoc solution in one dimensional sim- 
ulations that successfully matches observation has been 
to invoke a transition from sub-sonic to supersonic flame 
propagat ion at some cri tical density, a "delayed deto- 
nation" (Khokhlov1 ll991I) . Phenomenologically, this al- 
lows the star to expand first before the detonation front 
runs through it, thus creating a wider region of incom- 
plete burning where Si, S, Ar, and Ca are produced be- 
fore the flame is quenched. Considerable effort is going 
into simulating the actual turbulent burning front over 
a range of different scales to look for evidence of such a 

4 This early quenching of the flame in one-dimensional defla- 
grations is characteristic of realistic flame v elocities that have su - 
perceded the "fast deflagration" model W7 of No moto et al.l (1984'l. 



transition llRooke et al.l2003URopke fe Hillebrandtl2004 
l2005HBell et al.ll2004albft . However, although the results 
of multidimensional physical effects and micro-physics of 
the flame can be worked into the one-dimensional models, 
the flame propagation is inherently three-dimensional. 
Current three dimensional simulations have generally re- 
frained from artificially inducing a transition to super- 
sonic burning, but these pure deflagrations still suffer 
from a deficit of intermediate mass products just as 
they did in one dimension. The other difficulty is that 
the three dimensional simulations are almost all well- 
mix ed, unlike the stratification seen i n Typ e la spec- 
tra l|Kozma et al.l 12005ft . iPlewa et all <H}04) have pro- 
duced preliminary models that retain a stratified ele- 
mental structure in three dimensions, by initiating the 
burning at a single, slightly off-center point and having 
the resulting hot bubble break out on the WD surface 
to trigger a 'gravitationally confined detonation'. How- 
ever if the seeds for instabilities are present at the earliest 
times they will rapidly run away t o thoroughly mix burnt 
and unburnt material (Rdpke & Hillebrandt 2005). 

The study of young SNRs is playing an increasingly 
important role in uncovering the SN explosion mech- 
anism. For example, the "jet" in the northeast cor- 
ner of Cassiopeia A (Cas A), first detected in the op- 
tical fast-moving S-rich knots by van den Borah & Dodd 
(1970), has been cited as a possible exa mple of the jet- 
induc ed core-coll apse explosion model ifKhokhlo v et alJ 
11999ft . However, iHughes et all l|2000aft found that the 
composition of the X-ray emitting portions of the jet was 
more easily explained by Si-rich explosive Oxygen burn- 
ing with little or no Fe, suggesting that the jet material 
originated further out in the exp l oding s tar. In contrast, 
to the southeast, Hug hes et al.l l)2000aft found knots of 
almost pure Fe out past the bul k of the Si-rich ejecta. In 
deeper Chandra observations, Hwang & Laming (2003) 
confirmed the composition of these knots, finding oth- 
ers that could be explained by a pure Fe plasma. The 
high-velocity Fe-rich bullets in Cas A are reminiscent of 
the bullets of near-core material that would be ejected 
as a convection cycle broke through the stalled shock in 
the ne utrino-driven explosion models of Bu rrows et alJ 
( 1995). This model was motivated by the early detection 
of Fe in the spectrum of SN 1987A. However, reproducing 
both features, Nc-rich bullets and a "jet" that originated 
further from the core, in one model is somewhat prob- 
lematic l|Burrows et al.ll2005ft. 

Type la SNRs have also shown interesting structures 
that lead to constraints on the explosion mechanism. 
Specifically, the stratification observed in Type la SN 
spectra is also seen in Type la SNRs. In the rem- 
nant of SN 1006 multiple lines of evidence point to 
stratification of the Si and Fe ejecta b oth from ultra- 
violet absorption lines of Si and Fe (e.g. [Hamilto n et alJ 
1997) and the l ack of Fe-L emission in the X-rays 
IjVink et all 12003ft . In Tycho's SNR, spatially resolved 
spectroscopy of the ejecta supports the idea of strat- 
ification of conditions and composition llVancur a et alJ 
ll995llHwang fc Gottheldll997llBadenes et alJl2003ft . De- 
tailed comparisons of Tycho to one and three dimen- 
sional Type la SN explosion models have shown that one- 
dimensional, stratified "delayed detonation" models best 
reproduce the spatial composition of Tycho's SNR while 
recent well-mixed three-dimensional deflagration models 




Fig. 1.— ATCA 1.4GHz (left) and 5GHz (right) images of SNR G337.2-0.7 smoothed to 15"resolution (North:up Eastdeft). 
intensities are linearly scaled from to 20 mjy per beam (1.4GHz) and to 10 mjy per beam (5GHz). 
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fail <|Badenes et alJl2005b|) . 

These considerations regarding SN explosions are 
the context for our interpretation of the new XMM- 
Newton, Chandra and ATCA X-ray and radio obser- 
vations of G337.2— 0.7, an ejecta-rich remnant that still 
bears a clear imprint from the explosion that produced it. 
In Section 2 we outline the radio and X-ray observations 
and data reduction. The radio and X-ray morphologies 
of G337.2— 0.7 are compared in Section 3 and an estimate 
of the distance to the remnant is found from H I absorp- 
tion measurements in Section 4. Section 5 discusses the 
results of non-equilibrium ionization (NEI) modeling of 
the global X-ray spectrum of the remnant. The valid- 
ity of utilizing a single temperature model for the global 
spectrum is investigated with respect to spatial varia- 
tions throughout the remnant in Section 6. We then 
proceed to compare the results of our global NEI model 
to the predicted abundances from core-collapse and Type 
la SN models in Section 7. A detailed comparison of the 
global spectrum to th e Type la model SNR spectra from 
iBadenes et all l|2003j) is performed in Section 7.1. We 
conclude by incorporating the results of these various in- 
vestigations into a single picture of the ejecta-rich SNR 
G337.2-0.7 (Section 8). 

2. DATA, PROCESSING, AND CALIBRATION 

2.1. Radio Observations 

We have obtained radio observations of G337.2— 0.7 
at multiple frequencies using the ATCA, a six-element 
apertu re synthesis telesco pe located near Narrabri, Aus- 
tralia ((Frater et al. 1992). Observations centered at fre- 
quencies of 4672, 4800 and 5312 MHz were carried out 
on 1999 July 19/20 in the 0.75D configuration, using a 
bandwidth of 128 MHz at each frequency and for a to- 
tal effective integration time of 6.5 hours. Subsequent 
observations were carried out on 2004 Jul 13 (in the 6A 
configuration) and 2004 Dec 18 (in the 1.5D configura- 
tion), each recording data simultaneously at 1384 MHz 
(using a 128 MHz bandwidth) and in the H I line (us- 



ing a 4 MHz bandwidth) for a total integration time of 
20.5 hours. To recover the full range of spatial scales 
in the H I line, we also made use of shorter interferom- 
eter spacings from the Southern Galactic Plane Survey 
(SGPS) l|McClure-Orifnths et al]l2005l) 5 . 

The ATCA data were edited and calibrated using 
standard techniques. Observations at 4672, 4800 and 
5312 MHz were then combined into a single image 
(hereafter referred to as the 5-GHz image), while the 
1384 MHz data were imaged separately (hereafter re- 
ferred to as the 1.4-GHz image). Both images were 
deconvolved using maximum entropy, smoothed with a 
Gaussian restoring beam, and then corrected for primary 
beam attenuation. Residual atmospheric gain variations 
were significant for the 5 GHz data; these were corrected 
using two iterations of phase self-calibration. 

The final 1.4-GHz image has a resolution of 8" x 6" 
and a sensitivity of 70 //Jy beam -1 , while the 5-GHz 
image has a resolution of 12" x 10" and a sensitivity of 
90 /iJy beam -1 . The resulting maps are shown in Fig- 
ure ^ and in both cases are expected to be sensitive to 
all spatial scales associated with the SNR. The corre- 
sponding flux densities of the SNR are 1.55 ± 0.05 Jy 
at 1.4 GHz and 0.93 ± 0.02 Jy at 5 GHz, in both cases 
using o bservations of PK S B1934-638 referenced to the 
scale of lRevnoldsl l)1994fl . This implies a spectral index 
a = —0.4 (Sj, oc v a ), which is well within the typical 
range for SNRs. 

For the H I observation, continuum emission was sub- 
tracted from the data in the visibility plane. The data 
were then imaged at a velocity resolution of 3.3 km s , 
using a spatial filter aimed at excluding H I emission on 
larger scales while retaining sensitivity to absorption on 
spatial scales corresponding to structure in the SNR. The 
emission cube was then deconvolved using the CLEAN 
algorithm, and restored to a resolution of 39" x 36". 

2.2. XMM-Newton and Chandra X-ray Observations 

° http://www.atnf .csiro.au/research/HI/sgps/queryForm.html 




Fig. 2. — X-ray and radio images of SNR G337.2— 0.7. Upper left: broadband 1 to 5 keV Chandra image, binned by four pixels, linearly 
scaled from 1 to 12 counts per bin. Upper right: combined XMM-Newton MOS1 and MOS2 image of G337.2— 0.7 from 1 to 5 keV. 
Square-root scaling was chosen in this image to emphasize the fainter diffuse emission, almost filling the extent of the radio remnant. Lower 
left: Full resolution 1.4 GHz ATCA image linearly scaled from to 4 mjy per beam. Lower right: Overlay of the broadband Chandra X-ray 
(red) and the 1.4 GHz ATCA radio (blue-green) images. An arbitrary scaling was chosen to show both faint (blue) and bright (green) radio 
emission against a smoothed Chandra X-ray image (Gaussian kernel width 3 four-pixel bins) shown in red. 



G337.2-0.7 was observed for 38.4 ks by XMM- 
Newton on 2001 Feb 18, observation ID 0087940101. 
Standard reduction and processing of the three EPIC 
datasets were carried out using available XMM-SAS 6.0.0 
software. Starting from the level 1 event files, the lat- 
est calibrations were applied with the "emchain" and 
"epchain" tasks. The events were then filtered to re- 
tain only the "patterns" and photon energies likely for 
X-ray events: patterns to 4 (single and double pixel 
events) and energies 0.2 to 15.0 keV for the PN, pat- 
terns to 12 (including triple and quadruple events) and 
energies 0.2 to 12.0 keV for the two MOS instruments. 
Events near hot pixels, outside the field of view or oth- 
erwise suspect are "flagged" within the event file and 
were eliminated using the most conservative screening. 
Within the good time intervals, no significant flares were 
detected, and the background count-rate was consistent 



with that measured in average quiescent times. The ef- 
fect of any undetected soft flares was minimized by the 
use of local background from the same observation that 
was extracted from an ellipsoidal annulus outside the su- 
pernova remnant. Spectra extracted from the PN, MOS1 
and MOS2 cameras were fitted simultaneously to max- 
imize statistical significance (after consistency between 
instruments was verified). Detector response files and 
effective areas for the extracted spectra were also gener- 
ated with the standard XMM-SAS 6.0.0 software. 

Chandra observed G337.2-0.7 for 48.8 ks with the 
back-side illuminated CCD ACIS S-3, beginning on 20 
August 2002, observation ID 2763. These data were re- 
processed using CIAO 3.1 and caldb 2.26. The events 
were first calibrated with the appropriate instrument re- 
sponse files using "acis_process_events," then filtered on 
grade, on energy, and to eliminate bad pixels and periods 



of high background. The time-dependent gain correction 
was applied using the contributed software "TGAIN" 
from Alexei Vikhlinin (this is the same correction that is 
now available as part of "acis_process_events" in CIAO 
3.2). All spectra were extracted with their own weighted 
response files using the "acisspec" task. Background con- 
tamination, primarily from the Galactic Plane region, 
was important for G337.2— 0.7. For simplicity, an an- 
nulus around the SNR was used to subtract local back- 
ground directly from each spectra (weighted by the ratio 
of the areas of the background and source regions) . Im- 
ages of G337.2— 0.7 from X MM- Newton and Chandra 
are shown in Figures [3 and El 

3. RADIO AND X-RAY MORPHOLOGY 

The 1.4 GHz and 5 GHz images of G337.2-0.7 are very 
similar when convolved to the same spatial scale (Figure 
0]). Both show faint diffuse emission from a 4.5' x 5.5' 
(diam eter) SNR. A wid er radio image of the Galactic 
plane l|Green et alJl!99flt) shows G337.2-0.7 to be in an 
isolated region off the plane. While not a diagnostic of 
the progenitor type, one would expect a larger number 
of SNRs generated by Type la SNe than by core col- 
lapse SNe in this region, because Type la progenitor sys- 
tems have more time to drift away from their birthplaces 
at low latitudes before the explosion than more massive 
stars do. The sharp edges that presumably demarcate 
the outer blast-wave appear to be somewhat squared off 
and only limb- brightened in a few regions. It is impor- 
tant to note that this is a genuine feature rather than an 
artifact due to the use of "clean boxes" — i.e., while the 
deconvolution process was constrained by the use of re- 
gion selection, in both cases the boundaries of the regions 
used lay well outside the sharp edges seen here. 

The most prominent feature at both radio frequencies 
is a bright elliptical ring 2.0' x 3.2' in extent in the south 
part of the SNR. Such a brightening in the radio must 
indicate an enhancement in the magnetic field, electron 
density or both. The ring is incomplete and clumpy, ap- 
pearing primarily as short bright arcs with one spoke into 
the center and a bright clump in the north that extends 
outside the main ellipse. However, with the possible ex- 
ception of the bright northern clump that may have a 
slightly flatter spectral index (~ —0.25), the ring does 
appear to be a coherent structure with a constant spec- 
tral index, —0.5 ± 0.1 and weak linear polarization at 
5 GHz (up to peaks of 10% fractional polarization at the 
ends of the major-axis). 

The Chandra and XMM-Newton (MOS1+MOS2) im- 
ages are shown in the top panel of Figure The higher 
resolution Chandra image has been scaled linearly to em- 
phasize the bright clumps of X-ray emission near the 
radio ring, while the deeper XMM-Newton image is 
square-root-scaled to show the faint diffuse emission that 
extends throughout most of the radio shell. The bottom 
panel of Figure [21 shows the 1.4 GHz image at full res- 
olution on the left and a 3-color overlay of the Chandra 
X-ray (red) and ATCA 1.4 GHz (blue+green) images on 
the right. In this overlay, the fainter radio emission is 
shown in blue while the bright radio ring has been em- 
phasized by displaying it in green (turquoise). Purple 
and white regions indicate an overlap between the X-ray 
emission and the faint or bright radio emission, respec- 
tively. 



Looking outside of the bright radio ring, a wide band of 
X-ray emission north-east of the ring lies where the ra- 
dio emission is weakest. A faint radio filament traces 
the brightest portion of this swath of X-ray emission 
(compare the 1.4 GHz image with the 3-color overlay). 
Elsewhere some amount of faint X-ray emission accom- 
panies most of the faint radio emission inside the rim, 
with the possible exceptions of the narrow region south 
of the bright ring and the western-most side of the rem- 
nant. Along the limb-brightened portions of the radio 
rim, there are suggestions of X-ray emission above the 
background with one clear identification of an X-ray fil- 
ament aligned with the radio rim in the south-east. This 
filament is brightest in the X-rays where it is faint in the 
radio, but they are clearly related since they share the 
same small indentation (located at the brightest part of 
the X-ray emission of the filament). 

As for the bright ring itself, there is one short arc of 
X-ray bright material aligned with the northeastern tip 
of the radio ellipse, beginning just as the radio surface 
brightness is dropping. In the southern portion of the 
ring, the arcs of bright radio emission extend further out 
than their X-ray counterparts, while in the north-west no 
X-ray emission is seen from the strong radio arcs. The 
radio clump in the north of the ring shows only faint X- 
ray emission. Vice- versa the X-ray bright clump on the 
north-eastern side aligns with a particularly faint spot 
in the 1.4 GHz image, and the spoke of X-ray emission 
across the minor axis of the ellipse does not align well 
with the radio protrusion from the south rim. 

The meaning of the differences and correspondences 
between the X-ray and radio morphologies is not obvious, 
but clearly deserves more attention. What insights into 
the origin of the bright radio ring could the X-ray spectral 
observations divulge? If the ring demarcates the current 
position of the reverse shock, that would imply higher 
abundances in the X-ray spectrum from that region. An 
interaction with some structure in the circumstellar or 
interstellar medium (ISM) should be reflected in lower 
abundances of high-Z elements. Some feature intrinsic 
to the explosion itself would be likely to have different 
abundance ratios of heavy elements than elsewhere. We 
address these issues further in Sections 6 and 7.1. 

4. CONSTRAINTS ON THE DISTANCE TO SNR G337.2-0.7 

To constrain the distance to SNR G337.2-0.7 via H I 
absorption, we examined both the ATCA H I cube cen- 
tered on G337.2— 0.7, plus SGPS H I data on surround- 
ing fields, both on a channel- by-channel basis. Real ab- 
sorption features to G337.2— 0.7 should be significantly 
deep and conform to the shape of the remnant. Fig- 
ure OH plots the fractional absorption for three targets in 
the field, G337. 1—0.2 (an unresolved bright H II region 
~ 30' northwest o f the SNR, at a distance of 11 kpc 
ICorbel et al.l l|1999l) L the supernova remnant itself and 
a blank-sky position about 2' southwest of the remnant. 
The strongest absorption features toward G337. 1—0.2 are 
at Local Standard of Rest (LSR) velo cities of —22 km s - 1 
and —116 km s _1 (see also Fig. 4 of lSarma et alj fl997h 
The channel images show both to be entirely point- 
like and thus truly along the line of sight to this un- 
resolved H II region. SNR G337.2— 0.7 shows absorp- 
tion at —22 km s _1 which conforms to the shape of the 
remnant, but the next strongest absorption feature at 




Fig. 3. — H I absorption spectra from three regions in the field, 
G337. 1—0.2 an unresolved H II region ~11 kpc distant, the super- 
nova remnant itself, and a blank-sky position about 2' southwest 
at RA 16:39:14.151, Dec -47:53:43.52 (J2000). All absorption fea- 
tures to the H II region are point-like and hence real. Only the 
—22km s _1 feature in the remnant spectrum distinctly followed the 
morphology of the SNR. Absorption near —113km s — x extended 
in a wide swath a few arcminutes in size. The bottom panel shows 
a sample spectrum from a source free region 2' southwest of the 
remnant exhibiting clear absorption at —110km s _1 . Fractional 
absorptions relative to the 1.4 GHz continuum are plotted. Note 
that for the SNR these exceed an absolute value of 1 indicating the 
presence of H I self-absorption over this region, such as might be 
expected from the large molecular cloud seen in 12 CO. 



~ —116 km s _1 extends evenly across a wide swath of 
the image as can be seen in the blank-sky spectrum in the 
bottom pane l of Figure^ Given that the 12 CO survey of 
Dam e et alJ (J2001) shows that cold absorbing gas at this 
velocity covers the sight-lines toward both G337.2— 0.7 
and G337. 1—0.2 we take the lack of absorption specific 
to the SNR at —116 km s _1 to be significant and to in- 
dicate that the SNR is closer than this cloud. 

We thus adopt —22 km s _1 and —116 km s _1 as 
upper and lower limits, respectively, on the systemic 
velocity of G337.2— 0.7. Using the Galactic rotation 
curve of iFich et alJ 1)19891) . and assuming a distance to 
the Galactic Center of 8.5 kpc, we constrain the dis- 
tance to G337.2— 0.7 to be between 2.0 ± 0.5 kpc and 
9.3 ± 0.3 kpc, where we have assumed typical uncertain- 
ties of ±7 km s" 1 in the velocity of each H I feature. 
In the final analysis we also consider the possibility that 



general absorption in the region around G337.2— 0.7 has 
obscured a true absorption feature at —116km s _1 , i.e. 
that we should ignore the upper limit on the distance. 

5. X-RAY SPECTRAL ANALYSIS: THE GLOBAL 
SPECTRUM OF THE REMNANT 

5.1. An Aside on the Usage of the Word Abundance 

Given that the primary subject of this paper is to rec- 
oncile the abundance pattern that we see in the spectra 
of SNR G337.2— 0.7 with its radio and X-ray morphol- 
ogy, as well as directly with SN explosion model predic- 
tions, it is worthwhile to first clarify exactly to what we 
are referring. When speaking of the abundances there 
are many things to consider. The first distinction to 
make is between "fitted abundances," i.e. how much of 
any given element is inferred from a spectrum, and the 
overall metal production predicted by a SN model. The 
fitted abundances generally assume a single electron tem- 
perature and timescale and are actually proportional to 
the "volume emission measure" (EM) which is the prod- 
uct of electron density times the ion density integrated 
over the emitting the volume. Not only might some of 
the ejecta be at lower temperatures than others but if 
the density in a particular layer is too low it simply will 
not be seen. For example, in the Type la SNR models 
of Badene s et alJ ((2003) there are cases where the rela- 
tive emission measures of elements compared to relative 
masses produced can be different by 3 or 4 orders of mag- 
nitude. In fact, the less abundant element can appear to 
be th e more abundant on e (see O versus Fe in Figure 
4b of iBadenes et al.ll2003j) . Under these circumstances, 
it is clear that comparing globally fitted abundances to 
the bulk chemical composition obtained from a SN ex- 
plosion model is not always justified, and can lead to 
erroneous conclusions. The second hurdle is the nota- 
tion used for fitted abundances. Conventionally in NEI 
modeling, especially in XSPEC, the "abundance" of an 
element is defined by a linear factor relative to its "solar 
abundance." The solar abundance set is written on the 
logarithmic "dex" scale for the number of atoms (nx) 
of each element (X) relative to the number of hydrogen 
atoms (bh) as follows: dex(X) — 12.0 + log(nx/n~R)- In 
contrast, SN modelers generally report the solar masses 
of each element produced, not the number. Once the 
presence of ejecta in the emitting region has been de- 
tected by the fitted abundances, it is often the case that 
the abundances of the metals relative to each other are 
better constrained than the number of atoms relative to 
hydrogen. 6 Thus, throughout this paper we will refer to 
both the "abundance" of any given element (relative to 
hydrogen relative to solar) and the "relative abundances" 
or "abundance ratios" of one metal to another. These rel- 
ative abundances will be given with respect to their solar 
ratio. For example, a Ca to Si ratio of 2 times solar, does 
not imply that there are twice as many Ca atoms as Si 
atoms but rather that there are twice as many Ca atoms 
per Si atom as would be present in a "solar abundance 
plasma." Further, we will refer to the set of metal abun- 
dance ratios as the "abundance pattern" (listed relative 
to Si). Normalizing elemental abundances to a given "so- 

6 The line strengths are measured relative to the bremsstrahlung 
continuum which is generated by free electrons that are contributed 
by all elements in the plasma 




TABLE 1 
NEI Parameters for the Global Spectra of SNR G337.2— 0.7 
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Fig. 4.— XMM-Newton EPIC PN, MOS 1 and Chandra ACIS 
S-3 spectra of the overall SNR G337.2-0.7. (The MOS 2 spectrum 
is similar to MOS 1). The model is a non-equilibrium ionization 
thermal plasma model with the abundances of Ne, Mg, Si, S, Ar, 
Ca and Ee allowed to vary along with the absorbing column den- 
sity, the ionization timescale and the electron temperature. An 
additional power-law component significantly improved the fit. An 
example of the contribution of the second component to the overall 
model is plotted here for MOS1. The spectra have been re-binned 
to signal-to-noise of 4 for plotting purposes only. 



lar" composition is useful for studies of the ISM, but it is 
often an awkward convention for SN ejecta, especially in 
the case of Type la SNe, which have no H in the ejecta. 

5.2. NEI Model Fitting Procedure 

The global SNR spectrum contained 13900, 14500, 
23900 and 34200 counts from MOS1, MOS2, PN, and 
ACIS-S3 respectively. To begin, the spectrum from each 
instrument was fit independently with a planar shock 
NEI model. The s hock model used is described in 
Hugh es et al.l (|2000bf) and allows a variable electron tem- 
perature, kT e , final ionization timescale, n ^t, and nor- 
malization. 7 This model includes lines from lMewe et al.l 
( 1985) that are impor tant during the t r ansien t ionizing 
stage as described in iHughes fc Singhl {1994). Planar 
shock models are widely used as a first approximation to 
the spectrum of a complex SNR in order to obtain esti- 
mates of the bulk properties and abundances of the ejecta 
(or blast-wave) by initially assuming that all emitting 
material is under the same conditions. During the fitting 
procedure, the abundances of all elements were initially 
held at their solar ratios, but freed as necessary where the 
data could clearly constrain the individual abundances. 
The NEI model was convolved with an absorption model 
(using solar abundances) for the interstellar column den- 
sity, N-r. Where the statistics warranted it, a power-law 
component was also included to account for additional 

7 The model was incorporated for our own use as a "local 
mode l" into XSPEC version 11.3.1 for spectral fitting (Arnaud 
1996, http://xspec.gsfc.nasa.gov/ I. 
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high-energy continuum emission that could be coming 
from either non-thermal X-ray synchrotron emission (ap- 
propriate for the power-law model used here) or a second, 
higher temperature thermal component (explored further 
in Sections 7 and 7.1). Although small shifts in the line 
centroids between instruments could be seen (of order 15 
keV, similar to those found in XMM-Newton analysis of 
Kelper's SNRlCassam-Chc nai et aLll2004h no significant 
differences in the fits between instruments were found 
for spectra extracted from any region. All small differ- 
ences in the fitted parameters could be explained by the 
different quantum efficiencies for the PN, the MOS, and 
ACIS-S3 which effected what portion of the spectrum 
had the greatest statistical weight. 

The PN, MOS1, MOS2 and ACIS-S3 spectra of the 
overall SNR were then fitted simultaneously with a sin- 
gle NEI model. The only parameter that was allowed 
to differ between the instruments was the normalization 
which varied slightly, most probably owing to imperfect 
area calculations around the chip gaps, hot and dead pix- 
els. An iterative process of fitting each of the instruments 
individually and then in combination was used for speed 
to locate the area of parameter space where the best fit 
would lie. This was followed by a thorough search of 
the surrounding area of parameter space using all four 
instruments to identify the best-fit and constrain the la 
uncertainty ranges of all the parameters. The results of 
these fits are given in Table Q] 

5.3. Results from the Global SNR Spectrum 

First let us consider how the overall SNR spec- 
trum compares with our previous ASCA results 
and conclusions abou t the remnant's dynamical state 
(Rakowski ct afl l2001l) . The overall best-fit column den- 
sity, electron temperature, ionization timescale and nor- 
malization derived from the ASCA spectrum are all 



reasonably consistent with the more well-constrained 
XMM-Newton and Chandra results. The main result 
from the ASCA spectrum was the presence of a signif- 
icant amount of ejecta as exhibited by the prominent 
Si, S and Ar lines. The overabundance of these three 
species is clearly demonstrated in the new observations, 
with a better and more tightly constrained fit to their 
ionization state, allowing us to more precisely determine 
their abundances. Furthermore, a strong Ca component, 
which was hinted at with one outlying point above the 
ASCA model, is firmly detected with the far superior ef- 
fective area of XMM-Newton and Chandra. The only 
difference between the ASCA and the XMM-Newton 
and Chandra best-fit abundances lies in the prediction 
for Ne, Mg and Fe. The ASCA results for these species 
were based solely on the f-2 keV portion of the spectrum, 
where Ne and Mg Ka lines and Fe-L lines should lie. 
However, no strong evidence for particular lines was seen, 
and the model preferred to fill in that region with Fe- 
ll/ emission. The XMM-Newton and Chandra spectra, 
however, exhibit a prominent Mg line and should be sen- 
sitive enough to reveal an Fe-K line (at ~6.7 keV), under 
the assumptions of a single temperature planar model, 
if Fe were truly super-solar in abundance. Hence the 
deeper XMM-Newton and Chandra spectra prefer the 
combination of a less-absorbed, cooler continuum plus Ne 
and Mg lines rather than a more highly absorbed contin- 
uum and Fe-L lines to fit the soft emission. However, 
even this best-fit model only approximates, but does not 
reproduce, the features around the He-like Mg K-shell 
"triplet" at ^1.34 keV. For any single instrument, dis- 
crepancies from the model of similar magnitude occur 
across the entire spectrum (see the residuals in Figure 
HJ. However, for the most part the residuals from differ- 
ent instruments do not all agree in sign or trend, whereas 
around 1.5 keV and below all three exhibit a similar pat- 
tern (the MOS 2 spectrum is similar to that of MOS 1). 
Combined with the fact that there are two alternative 
minima in \ 2 to explain the soft emission, it is fair to 
say that the Ne, Mg and Fe results are less conclusive 
than for Si, S, Ar and Ca. 

The new X-ray observations strongly indicate the pres- 
ence of an additional spectral component. The greater 
than 150 reduction in x 2 for the addition of a powerlaw 
is highly significant (F-statistic probability of chance oc- 
curence: 5.6xl0 -20 ). Unfortunately the nature of the 
second component is not constrained by these spectra. 
An additional thermal component results in an equally 
viable fit (see Section 7). For the purposes of this sec- 
tion the question is how has the addition of a powerlaw 
component changed the fitted abundance pattern from 
the XMM-Newton and Chandra spectra. As would be 
expected, the abundances increase systematically when 
a second component is allowed to explain part of the 
continuum emission. However the relative abundances 
of Si,S, and Ar are virtually unchanged and [Ca/Si] is 
increased but only from 3.09 to 3.4 [Ca/Si]©. Without a 
second component, the Ne and Mg abundances relative 
to Si drop below solar, but their presence is still required 
for the fit. The Fe abundance in the single component 
fit drops steeply partially because a higher electron tem- 
perature (1.1 versus 0.74 keV) is required to explain the 
high energy continuum. 

Given our new X-ray parameters along with the more 



well-defined size and distance from the radio observa- 
tions we can refine our previous estimates of the ambient 
density and age of the remnant to see if these are con- 
sistent with a young ejecta-dominated SNR. For a first 
estimate, we assume that the normalization and tem- 
perature of our single temperature fits are dominated 
by the blast-wave component so that we can compare 
them to an analytical Sedov solution for the expansion 
of a middle-aged SNR into a uniform density medium. 
Following Hughes ct al. (1998), we obtain the following 
values for the ambient density, nn, age, explosion en- 
ergy, Eq, and amount of swept-up interstellar material, 
Msuy given the electron temperature, kT c , normaliza- 
tion N = n e nnV/4irD 2 , maximum angular radius Or, 
and an upper limit on the distance to the remnant D: 
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The upper limit on the distance implies that the rem- 
nant has reached the Sedov stage, and that the spectrum 
should be dominated by the outer-blast- wave moving into 
the ambient medium. On the other hand, if the lower 
limit on the distance is used, the remnant is still quite 
young, ~ 750 years, and the blast-wave would only have 
swept up O.8M , but suggests an implausibly low ex- 
plosion energy, Eq = 1.0 x 10 49 ergs (even younger and 
less energetic for the minimum angular radius). There 
are two other potential indicators of the age of the rem- 
nant, the abundances relative to solar and the ionization 
timescale. High-Z elements such as Ca have been shocked 
and yet the absolute abundances are less than a factor of 
ten. The ionization timescale measured is also longer 
than predicted by either end of the allowed distance 
range. Both of these would seem to prefer older ages 
with a greater amount of swept-up material. However 
they could also be explained by the reverse shock hav- 
ing only reached a small portion of (well-mixed) ejecta, 
which if it were hot would bias the single temperature fit 
toward a longer ionization timescale. Further discussion 
of the evolutionary state of G337.2— 0.7 and the inter- 
play of ejecta and ISM contributions to is deferred to 
Section 7.1 when interpreting the detailed Type la SNR 
models that explicitly incorporate the ejecta contribu- 
tion. 

Our NEI fits require that the rat io of Ca to Si is at 
least 3.4__q;75 times the solar value (An ders fc Grevessa 
1989) at the 99% confidence level. However, this is only 
true if our model assumptions hold; i.e. that the SNR 
can be adequately described by a single kT c , single fi- 
nal n e t, planar shock. Density and composition profiles 
in the ejecta, instabilities both in the explosion and in 
the SNR evolution, and profiles or inhomogeneities in the 
ambient medium, should translate into different temper- 
atures, densities, and timescales as a function of posi- 
tion (and composition) throughout the remnant. Such 
variations are suggested by the clumpy morphology of 
G337.2— 0.7 in the X-ray band, an anchor-shaped bright 




Fig. 5.— Left: Three-color Chandra image of SNR G337.2-0.7: red, 1.0 to 1.68 keV, green, 1.68 to 2.3, blue, 2.3 to 5.0 keV (binned by 
4 pixels, linearly scaled from zero to half the maximum number of counts per bin in that band). These bands were chosen because they 
showed the largest variations from each other. The extraction regions chosen for spectral comparisons are as marked. Right: A three-color 
image from XMM-Newton MOS1 and MOS2, using the same three energy bands, but square-root scaling. 



region to the south surrounded by smoother, more dif- 
fuse emission. To test the adequacy of our single planar 
shock model, we first investigated how far individual re- 
gions deviate from the globally averaged spectrum. 



6. 



X-RAY SPECTRAL ANALYSIS: SPATIAL VARIATIONS 
WITHIN THE REMNANT 



A comparison of images in small energy bands was used 
to identify possible areas of spectral variation. Spectra 
from these regions were first modeled independently of 
the global spectrum, each with a single temperature, sin- 
gle timescale, variable abundance planar shock model. 
Even in these sub-sections of the remnant, the planar 
shock models might not be representative of the real 
conditions in the ejecta, but the uniformity of model as- 
sumptions allows for fair internal comparisons between 
regions and to the global fitted parameters. The indi- 
vidual variable-abundance fits are used to compare the 
electron temperatures, ionization timescales and [Ca/Si] 
ratios throughout the remnant independent of the global 
fit. Alternatively, we can test if plasma conditions alone 
are sufficient to explain the variations by fixing the rela- 
tive abundances to the values found in the global fit and 
allowing only the temperature, absorbing column den- 
sity and ionization timescale to vary. The aim of this 
search is twofold, first to identify inhomogeneities in the 
SNR, but also to test if modeling the sum of these dis- 
parate regions with a single temperature, timescale and 
abundance set has biased our results. 



6.1. Extraction Regions 

extracted from the entire low-surface- 
region show no obvious spectral diffcr- 



Spectra 
brightness 

ences from a spectrum containing all the high-surface- 
brightness clumps, hence surface-brightness by itself is 
not a good criteria for identifying spectrally distinct re- 
gions in G337.2— 0.7. Instead, narrow energy-band im- 
ages from each instrument were examined separately to 
identify any spectral variations across the SNR. Each 



event-file was divided into small spectral bins covering in- 
dividual lines and areas of continuum emission identified 
in the overall SNR spectrum. Images in these bins were 
compared to look for energy bands whose images differed 
strongly from each other. Energy bands that were indis- 
tinguishable spatially were combined for greater signal to 
noise. The three energy bands that showed the greatest 
amount of variability from each other are displayed in 
Figure E 1.0 to 1.68 keV (red) 1.68 to 2.3 keV (green) 
and 2.3 to 5.0 keV (blue). This three-color image of 
G337.2— 0.7 separates the energy range rich in Si and 
S lines (green) from the harder, more continuum domi- 
nated, regime (blue) and the soft band (red) which may 
include contributions from Ne, Mg, or Fe-L. The largest 
fluctuations in color are within the bright regions, but 
similar variability in the more diffuse region should not 
be ruled out, since it might simply be missed at this level 
of signal to noise. The differences seen within the bright 
region must be real features since they are repeated in 
both the XMM-Newton and Chandra observations. 

Eleven regions were chosen on the basis of their 
colors in Figure [SJ "bluel," "blue2," "greenl," 
"green2," "riml," "rim2," "whitel," "white2," "white3," 
"yellowl," and "yellow2." Spectra were extracted from 
these regions for each instrument, and individual re- 
sponse matrices and ancillary response files were con- 
structed using standard XMM-SAS and CIAO tools. 
Their spectra are shown in FigureEl All exhibit strong Si 
and S lines indicative of enhanced abundances, but vari- 
ations in conditions are evident in the relative strength 
of the He-like line complex and the H-like line of each 
species. For instance "green2" shows little or no evi- 
dence of any H-like Si, as opposed to "whitel" where the 
2.006 keV line from H-like Si is relatively strong. Because 
the regions were chosen to isolate potential spectral dif- 
ferences indicated by our image analysis, they do not all 
have similar numbers of counts. Hence not all spectra 
can constrain the abundances equally well. 
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Fig. 6. — Spectra of eleven small regions around the SNR. In the NEI models shown here, the elemental abundance ratios were fixed to 
those found in the global SNR fit, without an additional power-law component. The spectra from all three XMM-Newton instruments 
and Chandra ACIS-S3 were used in constraining the best fit model, but only the XMM-Newton EPIC PN spectra, which have the highest 
signal-to-noise, are plotted here for clarity. 



As before, the spectra were fit with a single component 
NEI model, with a variable electron temperature, ioniza- 
tion timescale and column density. For each region we 
first determined which elements required non-solar abun- 
dances as measured by the F-test for the addition of a 



new parameter. All regions showed significant signs of 
Si and S enhancement, and many (bluel, blue2, greenl, 
green2, whitel, white2) required an overabundance of 
Ca. No regions had statistics sufficient to determine the 
Ne, Mg, Ar, or Fe abundances compared to their solar 
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values or constrain an additional power-law component. 
However, some do appear to have stronger Mg lines than 
others. 

6.2. Two Approaches and their Results 

In order to directly compare abundance measurements 
between regions, their spectra must be modeled under 
identical assumptions. Since the most interesting results 
from the global fits were those of Si, S and Ca, we chose 
to fit all regions with only those abundances free (all pa- 
rameters, Table |21 selected values repeated in Table OJ. 
Overall, the fitted abundances of Si, S, and Ca are super- 
solar but tend to be slightly lower than in the global 
spectrum — probably because no power-law component 
was included in the small region fits. At la the fitted 
parameters indicate variation in the electron tempera- 
tures (0.61+E5;£7 keV "riml" to 1.17 ± 0.04 keV "bluel"), 
logarithmic timescales (11.09io"i5 log(cm~ 3 s) "rim2" to 
fully equilibrated, 13.0_o 8 log(cm _3 s) "white2"), and 
the [Ca/Si] abundance ratio (1.9±0.5 "bluel" to 4.0±0.8 
"white2"). The detection of variations in the plasma pa- 
rameters holds at a much higher significance level than 
that of the abundances. Also note that the fitted temper- 
atures and [Ca/Si] ratios do not appear to be correlated 
given that 7 out of 11 best-fit [Ca/Si] ratios lie between 
3.69 and 4.22 over a wide range of temperatures. From 
these fits it is clear that inhomogeneities in the conditions 
of the emitting material exist. 

The question at hand is whether the highly super-solar 
[Ca/Si] ratio found in the best-fit to the global spectrum 
could somehow be an artifact of naively summing spectra 
from regions with varying plasma conditions. The spec- 
trum from region "bluel" indicates that there do exist 
areas in the remnant where the [Ca/Si] ratio is consis- 
tent with the solar ratio and inconsistent (at la) with the 
high [Ca/Si] ratio found from the global spectrum. Other 
than this one region, the spectra are consistent with the 
conclusions from the global model but are too shallow to 
confirm it with any degree of certainty. At best, three re- 
gions, "whitel," "white2" and "green2" have [Ca/Si] ra- 
tios that exceed 3 times solar at the la level. Hence these 
fits have only shown that high abundances of Ca relative 
to Si may exist throughout most of the remnant, but have 
not substantiated the argument. On the other hand, nei- 
ther did we find a correlation between high [Ca/Si] ratios 
and high temperatures or a single region especially rich 
in Ca compared to Si that could account for the high 
Ca-K line flux in the overall spectrum. 

Alternatively we can model each of the individual re- 
gions with the relative abundances found for the overall 
spectrum to see if any of the spectra require abundances 
significantly different than the overall fit or if all the vari- 
ations can be fully accounted for by changes in the con- 
ditions. If we fix the relative abundances of Ne, Mg, Si, 
S, Ca, Ar, Fe, and Ni to the best-fit values for an NEI 
model of the global spectrum (without a power-law com- 
ponent) but allow the overall metallicity to be free, the 
spectra from all but one region ("yellow2") are as well 
fit as they were under the previous set of assumptions 
(Table |2J . Examining the spectrum of the discrepant 
region, "yellow2," in Figure it appears that much of 
the disagreement with the overall SNR ratios is at Ne 
and Mg. Using the F-statistic for the addition of 2 free 



parameters between the model which assumes the over- 
all SNR abundance pattern and that which allows Si, S, 
and Ca to vary freely, gives a 99.98% probability that 
the allowing for independent abundances of Si, S, and 
Ca is a true improvement for this region. Given that the 
fixed abundance set used is not independent of the small 
region spectra, the F-statistic is not strictly correct and 
actually underestimates the significance of the reduction 
in x 2 . 

Given that the two prescriptions above identified two 
different regions as discrepant, we see that the conclu- 
sions of any given analysis are effected by the choice 
of how to treat all the abundances, not just the abun- 
dances of the elements with the most prominent lines. 
The "independent" fits which allowed Si, S, and Ca to 
vary but held all others at their solar values showed that 
the [Ca/Si] ratio for "yellow2" is entirely consistent with 
the overall SNR result despite the poorer fit around the 
Ne and Mg lines when one requires all abundances to fol- 
low the SNR pattern (as compared to the freely varying 
Si, S and Ca model). Vice versa, while the "independent" 
fits indicated that "bluel" has a lower [Ca/Si] ratio than 
the overall SNR, this spectrum is marginally better fit, 
Ax 2 = 8.0, assuming the overall abundance pattern of 
the SNR rather than forcing some elements to be fixed 
to their solar abundances. In the case of "yellow2" the 
abundances do not match the overall SNR, however it is 
not the high [Ca/Si] ratio that is the problem. In the 
case of "bluel" the goodness of fit to the overall SNR 
abundance ratio model shows that the initial result of 
a lower [Ca/Si] ratio was an artifact of not properly ac- 
counting for the emission from the less prominent species. 
However, it is not just the abundances themselves that 
are affected by the starting assumptions. The tempera- 
tures found using the global SNR abundance ratios are 
systematically higher by about 0.2 keV than when only 
Si, S, and Ca are freed from their solar values. This 
demonstrates the magnitude of inherent bias that any 
given assumption about the abundances can make in an- 
alyzing spectra with insufficient counts to constrain all 
the relevant abundances. The solution to this problem 
of a strong dependence on model assumptions would be 
a full Bayesian analysis, but G337.2— 0.7 would not be a 
good test case for such a technique. For the present, the 
reader is urged to treat abundance measurements with 
caution and to test the sensitivity of their results to their 
underlying assumptions. 

Despite variations in conditions, the average of the 
metallicities of the small regions compared to the overall 
SNR model is 1.01, or 1.14 when weighted by the total 
number of bins in each spectrum. Hence, at least in this 
instance, using a single global model has, if anything, un- 
derestimated the total abundances. Likewise the Ca to Si 
ratio averaged over the small regions is 3.96 times the so- 
lar ratio (or 3.77 weighted by the number of bins) which 
is also greater than the value of 3.09 times solar found 
in the global SNR fit before the inclusion of a power-law 
component. 

6.3. Implications for the Morphology of G337.2-0.7 

Disregarding for the moment the low significance level 
of any abundance variations, what would the fits to the 
small regions in Tables [21 and imply about the X-ray 
and radio features discussed in Section 3? The X-ray 



12 









TABLE 2 














Fits for small 


REGIONS IN G337.2 


-0.7 






Region 


(10 22 cm- 2 ) 


kT B 
(keV) 


log(n c t) a 
(log(cm — 3 s) 


Si 


S 


Ca 


[Ca/Si] 


riml b 


3.65 (3.41-3.94) 


0.61 (0.54-0.72) 


11.67 (11.11-12.32) 


1.50 (1.24-1.84) 


2.34 (1.93-3.03) 


11.76 (0-40.7) 


7.84 (< 30) 


rim2 b 


3.35 (3.21-3.55) 


0.82 (0.72-0.87) 


11.09 (10.94-11.26) 


2.19 (1.95-2.40) 


3.27 (2.89-3.86) 


9.25 (1.94-16.93) 


4.22 (0.9-8.6) 


white 1 


3.01 (2.90-3.09) 


0.93 (0.91-0.95) 


12.0 (11.77-12.56) 


3.04 (2.80-3.32) 


4.13 (3.80-4.55) 


11.2 (9.02-13.58) 


3.69 (3.0-4.5) 


white2 


3.36 (3.20-3.46) 


0.87 (0.85-1.02) 


13.0 (12.20-13.00) 


2.33 (2.13-2.68) 


3.61 (3.31-4.08) 


9.28 (7.34-11.30) 


3.98 (3.15-4.8) 


white3 b 


3.02 (2.92-3.17) 


0.90 (0.83-0.90) 


12.50 (12.17-13.00) 


2.10 (1.84-2.33) 


2.59 (2.29-2.87) 


5.86 (3.83-8.50) 


2.79 (1.8-4.2) 


bluel 


3.84 (3.69-4.03) 


1.17 (1.13-1.21) 


12.13 (11.98-12.46) 


2.78 (2.52-3.08) 


3.48 (3.13-3.86) 


5.28 (3.91-6.65) 


1.90 (1.4-2.4) 


blue2 


3.59 (3.37-3.76) 


1.08 (0.96-1.23) 


11.41 (11.25-11.57) 


2.53 (2.28-2.89) 


3.24 (2.94-3.75 ) 


10.05 (7.6-13.1) 


3.96 (2.9-5.2) 


yellow l b 


3.37 (3.17-3.54) 


0.67 (0.61-0.76) 


11.48 (11.17-11.80) 


1.32 (1.17-1.50) 


1.31 (1.10-1.55) 


5.40 (0-15.15) 


4.09 (< 12) 


yellow2 b 


2.79 (2.69-2.89) 


0.87 (0.81-0.95) 


11.49 (11.36-11.64) 


1.50 (1.37-1.65) 


1.47 (1.30-1.65) 


5.89 (2.61-9.38) 


3.93 (1.8-6.3) 


greenl 


3.24 (3.20-3.39) 


0.81 (0.74-0.85) 


11.51 (11.42-11.65) 


3.19 (2.93-3.39) 


4.60 (4.32-4.89) 


10.12 (7.76-12.69) 


3.17 (2.4-4.0) 


green2 


3.87 (3.68-4.03) 


0.77 (0.72-0.85) 


11.49 (11.34-11.60) 


3.63 (3.36-3.95) 


4.64 (4.25-5.11) 


14.60 (10.95-18.61) 


4.02 (3.0-5.2) 



a By log(n c t) = 13.0 ionization equilibrium has been reached, hence no upper limit can be made on the age or density for "white2" or "whitc3.' 
b the F statistic did not support the Ca abundance as an additional free parameter in these regions 



TABLE 3 

Small Regions modeled with overall abundance ratios 







solar abundances 






abundance ratios 






except 


variable Si,S,C 


!a a 


fixed at global SNR values 


region 


kT c (keV) 


[Ca/Si]/[Ca/Si] b 


X 2 (d.o.f.) 


kT c (keV) 


Y/Ysnr c X 2 (d-o.f.) 


riml 


0.61 




7.84 


121.11 (104) 


0.84 


0.65 123.61 (106) 


rim2 


0.82 




4.22 


308.10 (271) 


1.01 


1.02 307.42 (273) 


white 1 


0.93 




3.69 


305.18 (241) 


1.14 


1.29 309.59 (243) 


white2 


0.87 




3.98 


299.26 (212) 


1.16 


1.09 295.76 (214) 


white3 


0.90 




2.79 


167.92 (173) 


1.20 


0.95 165.65 (175) 


bluel 


1.17 




1.90 


220.03 (200) 


1.41 


1.12 212.05 (202) 


blue2 


1.08 




3.96 


290.07 (275) 


1.19 


0.90 289.53 (277) 


yellow 1 


0.67 




4.09 


100.04 (116) 


0.81 


0.42 100.12 (118) 


yellow2 


0.87 




3.93 


239.91 (205) d 


0.98 


0.43 261.75 (207) 


greenl 


0.81 




3.17 


530.43 (436) 


0.92 


1.55 531.24 (438) 


green2 


0.77 




4.02 


405.68 (398) 


0.844 


1.72 400.01 (400) 



a Selected values reproduced from Table |21 

b c.f. [Ca/Si]/ [Ca/Si] for SNR overall without an additional power law component is 3.09. 

c Ne, Mg, Si, S, Ar, Ca, and Fe were tied at their relative abundances from the global SNR fit. Y/YgjvK 
represents a ratio of the form [Nc/H]/[Nc/H]sjv.r for each of these species. 

d The F-statistic for additional parameters from the fixed abundance ratio model to the variable Si, S 
and Ca model is 9.33, probability of chance occurrence: 1.32xl0 — 4 . 



regions that most closely follow the bright radio ring in 
both position and surface brightness are "yellowl" and 
"yellow2." These two regions have the lowest abundances 
of Si, S, and Ca found in the remnant and also the lowest 
metallicity in the fits where the abundance ratios were 
fixed to the overall SNR values. Further, keeping Ne, 
Mg and Fe at their solar abundances produced a signif- 
icantly better fit for region "yellow2" than the sub-solar 
overall [Mg/Si] values. Hence it is plausible to conclude 
that there is more swept-up material in "yellowl" and 
"yellow2" than elsewhere and that the bright radio ring 
denotes an interaction with a denser feature in the cir- 
cumstellar or interstellar medium (such as a circumstel- 
lar wind if the remnant is young or an interaction with 
a molecular cloud if it is older) not the current position 
of the reverse shock. While the abundances are higher 
for regions "whitel," "white2," and "white3," the nearly 
equilibrated ionization timcscalcs support the hypothesis 
of higher density material around the ring. 

The results from the two "rim" regions are less con- 
clusive. While "riml" does have the third lowest abun- 
dances, it still favors a high Ca to Si ratio indicative of 



newly synthesized material. The larger region "rim2" 
shows no sign of lower abundances compared to the rest 
of the remnant, and the highest absolute abundances 
are found in the unremarkable, diffuse emission regions 
"greenl" and "green2." This may be indicating that the 
reverse shocked ejecta and the blast-wave into the ISM 
are no longer distinct in this remnant quite unlike other 
middle-aged SNRs such as DEM L71. On the other hand, 
if G337.2— 0.7 lies at the upper end of our distance limits 
it would be far from the solar neighborhood, closer to the 
Galactic center where considerably more recycled mate- 
rial should enhance the metallicity by as much as 0.3 dex 
(e.g. iShaver et~aIIIT983l iRolleston et aLW oOO). The sec- 
ond solution may seem simpler, but it exacerbates the 
problem of a super-solar Ca to Si ratio by requiring this 
abundance ratio to have been obtained over the multiple 
SN explosions that produced the metals in the current 
Galactic-center ISM. 

6.4. Summary 

In summary, the analysis of small regions through- 
out the remnant agree with the conclusions of metal- 
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enhanced material and a super-solar Ca to Si ratio from 
the overall spectrum. While temperature and ionization 
timescale differences are present, abundance variations 
are comparatively small, as measured by planar shock 
models. In particular, we do not find any evidence for a 
pocket of hot Ca that could explain the Ca-line flux in 
the overall spectrum, nor do we see any strong deviations 
from the overall SNR abundances except in a portion of 
the bright ring which appears to have a greater contri- 
bution from swept up material. 

7. POTENTIAL SN MODELS 

Core-collapse and Type la explosions have widely dif- 
ferent patterns of abundances that should be readily de- 
tected in the X-ray emission from their ejecta. In core- 
collapse explosions the explosive nuclear burning front 
passes through the existing hydrostatic layers of the star, 
burning them incompletely and leaving large quantities 
of intermediate products like O, Ne, Mg, Si, S, Ar and 
Ca. Much of the Fe-group material produced may fall 
back onto the compact object since these are formed 
at the highest temperatures and thus nearest the cen- 
ter, although some ca n be flung out to larg e distances as 
seen in Cassiopeia A IjHughes et al.ll200r)a^ . In Type la 
explosions, on the other hand, the thermonuclear burn- 
ing efficiently incinerates most of the C/O white dwarf 
to Fe-group elements, with fewer intermediate products 
and some unburnt C and O. There is no compact object 
left behind in Type la's, so all the Fe-group material is 
ejected. Hence, traditionally the O/Fe ratio has been 
used to identify SNRs as core-collapse or Type la ex- 
plosions (although the vastly different conditions of the 
O and Fe layer s in a Type la can m ake this fitted ratio 
deceiving; see Badcncs et al.ll2003l) . Here the question 
is whether a highly non-solar [Ca/Si] ratio could also be 
diagnostic of the explosion mechanism. 

In core-collapse SNe, Si, S, Ar and Ca are all pri- 
marily produced in incomplete explosive oxygen burning, 
where their relative abundances are set to near solar val- 
ues by quasi-equilibrium a nd do not vary sig nificantly 
with conditions or mass (Wo oslev. Arnett. fc Clavton 
1973). Highly energetic or highly asymmetric explo- 
sions have larger volumes of "a-rich freeze-out" which 
produces 44 Ti and 44 Ca. However, these explosions 
also have extended volumes of explosive oxygen burn- 
ing such that the common Si and Ca is otopes domi- 
nate the [Ca/Si] ab undance ratio (e.g. iMaeda et alJ 
2002; Kifonidi s et al. | |2003ft. For the p r ototypical core - 
collapse models of IWooslev fe Weaver! 11 9951 ( WW95) , 
iRauscher et"all 1J2002J) . and iLimongi fe Chicffi (2003), 
the highest [Ca/Si] ratio found (~F8|Ca/Si]©) was 
for the 12 M Q "S12A" model of 1WW951 hereafter 
WWS12A. 8 

For Type la models (e.g. iHoeflich fc Khokhlovl 1199ft 
ITwamoto et al.llT999l IBadenes et al .1120031) . high [Ca/Si] 
ratios are found in delayed detonation explosions (cur- 
rcntl y favored by the SN light-curves and spectra). 
In Badcn es et alJ l)2003|) . the [Ca/Si] ratios of a de- 
layed detonation, DDTe, and a pulsed delayed deto- 
nation, PDDe, model are 2.5, and 2.9 times the so- 

8 The total ejected mass (1.32.Mg) and kinetic energy of the 
explosion (1.17 X 10 51 ergs) of WWS12A are quite similar to the 
values for Type la explosions, making a comparison on this basis 
largely moot. 



TABLE 4 
Relative Abundances: [X/Si]/[X/Si]© 

Ratio G337.2-0.7 WWS12A a PDDe b DDTe b 



Ne/Si 

Mg/Si 

S/Si 

Ar/Si 

Ca/Si c 

Fc/Si d 



1.19 
0.50 
1.44 
0.49 
3.40 
0.28 



0.13 


0.0010 


0.015 


0.13 


0.0017 


0.025 


1.46 


1.5 


1.4 


2.20 


0.68 


0.60 


1.8 


2.9 


2.5 


0.37 


0.89 


0.91 



"•derived production fa ctors after radioactive decay: 
IWooslev & Weaver! JT553, table 6A). 

'"Badcncs et al. (2003) Ne and Mg masses were not pre- 
sented but are shown here. 

c 2.65 - 4.2 at 99% confidence. 

d 0.14 - 0.46 at 99% confidence. 



lar ratio, respectively. Given the stratification of the 
ejecta present in the one-dimensional models, the lack 
of observed Fe in G337.2— 0.7 can be explained by the 
fact that Fe is shocked last and stays at a lower den- 
sity (see iBadenes et al.l 120031 for a discussion) . Re- 
cent 3-D deflagration models of Type la explosions pro- 
duce a similar abundance pattern to the delayed deto- 
nation models shown here, without an a rbitrary transi- 
tion from deflagration to detonation fe.g.lReinecke et al] 
I2001 ICamezoet all 120(131 Iflarcia-Senz fc Bravol 12005). 
However, they su ffer from a sho rtage of intermediate- 
mass elements (IBravo et all 2005) and the burnt and un- 



burnt material are thoroughly bulk-mixed such that it 
would be difficult to explain the lack of Fe emission in 
the data. 

Table 21 compares the fitted abundance ratios to 
the total metal production from a core-collapse model, 
WWS12A and two Type la models, PDDe and DDTe (re- 
call the caveats to such a comparison from Section [OJ. 
All models reproduce the S to Si ratio well. The core- 
collapse model matches the Fe to Si ratio well and does 
produce some Ne and Mg, but strongly over-predicts the 
production of Ar while still falling shy of the measured 
[Ca/Si] ratio. The [Ar/Si] and [Ca/Si] ratios are espe- 
cially problematic because high production factors of Ca 
isotopes were accompanied by high production factors 
of Ar in all o f the solar abundance massive star mod- 
els of WW95. Both Type la models clearly match the 
Ca to Si ratio best and do not conflict as badly con- 
cerning the production of Ar. On the other hand, Fe 
is strongly over-predicted and Ne and Mg are strongly 
under- predicted by the Type la models. To determine 
the degree to which these discrepancies can be seen in 
the spectrum we performed fits of the global spectrum 
with the abundance ratios of C, N, O, Ne, Mg, S, Ar, 
Ca, Fe and Ni (relative to Si) fixed to the SN explo- 
sion predictions. Figure \7\ compares NEI models for the 
spectrum that use the abundance ratios relative to Si 
of WWS12A or a de layed detonation (DDTe) model by 
IBadenes et all l)2003|) . The deviations from the data for 
WWS12A are both more prominent and harder to ex- 
plain away. While either model produces a good fit to the 
data if Ar (WWS12A) or Fe (DDTe) are allowed to vary, 
the stratification of Type la models provides a natural 
motivation for reducing the amount of Fe seen, whereas 
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Ar, Ca, and Si should be primarily co-spatial in the ejecta 
of a core-collapse explosion. Furthermore, swept-up so- 
lar composition material will add Ne and Mg needed by 
both models but should not flip the sense of the Ar to 
Si ratio for WWS12A. Tests of a two-component model, 
one with solar composition (ISM) the other with freely 
varying abundances (ejecta) bear this out. If the tem- 
peratures of the two components are set to reasonable 
values (such as the blast-wave and EM-averaged Si tem- 
peratures from Section 7.1) the fitted [Ne/Si], [Mg/Si] 
and [Fe/Si] abundance ratios in the "ejecta" component 
drop steeply, while the sub-solar [Ar/Si] ratio and super- 
solar [Ca/Si] and [S/Si] ratios remain at values similar 
to the one component fit. Hence a strict comparison of 
the global abundance ratios clearly favors an origin as a 
Type la SN. If true, G337.2-0.7 would be only the third 
potential la SNR with a sign ificantly disturbed morphol- 
ogy after Ke pler's SNR fe.g.lCassam-Chenai et alJl2004 
and N103B l|Hughes et al.lll995HLewis i et alJl200l) . both 
of which also have controversial identifications. 

A core-collapse scenario with solar-like abundance ra- 
tios between Si, S, Ar, and Ca could still be viable if the 
strong Ca-K line is the result of differing plasma con- 
ditions rather than abundance. Line fluxes in general 
are strongly dependent on the electron temperature and 
ionization, hence if the temperature of the Ca-rich ma- 
terial were sufficiently higher than in the Si-rich ejecta 
this could also account for their unusual line ratio. A 
global fit allowing a separate temperature and ioniza- 
tion time-scale for Ca but holding [Ca/Si] = [Ca/Si] per- 
forms equally well as our single temperature fit with a 
free Ca abundance, with almost identical best-fit \ 2 val- 
ues. The electron temperature associated with the Ca 
emission was 1.0 keV as compared to 0.75 keV for the 
Si and other elements. However, no spatially distinct re- 
gion with a high temperature and high [Ca/Si] ratio was 
identified in the small region fits. On the other hand, this 
does not preclude the possibility of two unresolved or co- 
existing plasmas. For instance, one could consider the 
case of a two density medium in pressure balance, with 
a single elemental composition shocked at a single time. 
We looked at a model with two NEI components with the 
same abundances, their temperatures tied at a series of 
fixed ratios and the ionization timescales offset appropri- 
ately for pressure balance. This did result in a sub-solar 
[Ne/Si] ratio (~O.8[Ne/Si]0) and reduced the amount of 
Fe required by the fit (^0.05[Fe/Si]o) but did not change 
the pattern of Si, S, Ar and Ca abundances substantially 
(~1.5[S/Si] , ~O.5[Ar/Si] , ~3.O[Ca/Si] ). Hence, al- 
though this question is not well-defined, it seems unlikely 
that a multi-temperature plasma with an abundance pat- 
tern typical of a core-collapse explosion could explain 
the global spectrum of G337.2— 0.7 without considerable 
compositional differences from place to place. 

In the above discussion we have compared the global 
spectrum of G337.2— 0.7 to the total metal output of the 
SN explosion models with the addition of a few simple 
two-temperatures tests and concluded that the products 
of a Type la explosion are more consistent with the spec- 
trum. In reality the blast-wave into the ISM may con- 
tribute substantially to the total spectrum and the ejecta 
itself will not be a homogeneous mix of the nucleosynthe- 
sis products. Layering of the composition or large scale 
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Fig. 7. — Total ejected elemental masses compared to the global 
SNR G337.2-0.7 spectrum. Only the MOS 1 spectrum is shown 
for simplicity, although all four instruments were considered in 
the analysis. The original free-abundances NEI model is shown 
in black. For comparison the relative abundances were fixed to the 
appropriate ratios given the total ejected mass of each element in a 
core-collapse (WW S12A: blue) and a Type la (DDTe: red) model. 



anisotropies imply that not all products will be visible in 
the X-rays at all times, nor will they be under uniform 
conditions. A true comparison of the global SNR spec- 
trum to explosion models requires that the evolution of 
ejecta and shock fronts be followed from the time of the 
explosion until the present age of the remnant to pre- 
dict the X-ray spectrum. While no such model is avail- 
able for the variety o f core-collapse models, the models of 
Badcncs et al. (2003) provide the opportunity for exactly 
such a comparison for the case of a Type la explosion. 

7.1. Investigation of Possible Type la Explosion 
Mechanisms 

Many groups have produced grids of Type la 
mod els to compare to SN spectra and ligh tcurves 
fe.g. iHoeflich fc Khokhlovl Il99a llwamoto et all fl999^ . 

Badcnes et al. (2003) have expanded upon this work to 
enabl e compar isons to the remnants of Type la SNe. 
iBadenes et all lj2003|) followed the plasma conditions and 
ionization state within the ejecta of Type la SNe from 
the explosion itself to 5000 years later, and calculated the 
expected X-ray emission of the ejecta and the blast-wave 
for various times in the evolution of the remnant. Their 
network of Type la explosion models attempted to cover 
all currently debated explosion mechanisms from pure 
detonation (super-sonic flame propagation) to pure de- 
flagration (sub-sonic flame propagation), as well as a sub- 
Chandrasekhar mass explosion. Their models showed a 
strong evolution of the contribution from each element 
to the X-ray spectrum. Indeed, for some models Fe did 
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not begin to appear until late times. 

The structure of the ejecta as prescribed by each SN la 
model will evolve as the forward shock expands out into 
the ambient medium and a reverse shock is generated 
that propagates back into the layers of ejecta. The X- 
ray emission from the SNR will have contributions from 
both the outer blast wave (reasonably approximated by 
a Sedov solution after on the order of 1000 years) and 
the reverse shock ( followed in detail in the models of 
iBadenes et al.ll2003j) . In this simple one dimensional pic- 
ture only the outermost layers will be shocked at the 
earliest times, and thus begin to ionize and emit in the 
X-rays. Elements that are only present in deeper layers 
will show up later in the spectrum and experience differ- 
ent temperatures and densities than the earlier layers. 

For comparison with the spectrum o f G337.2— 0.7, 
for ev ery explosion model in the grid of Ba denes et al.1 
(2003), the evolution of the SNR was computed in a grid 
of three ambient densities, p am (0.2, 1.0, 5.0 xl0~ 24 g 
cm -3 ). From each of these SNR models, spectra from 
the shocked ejecta were extracted at four different ages 
(500, 1000, 2000, and 5000 years). The synthetic spectra 
of Badenes et al. are completely characterized by four 
parameters only: the Type la explosion model, the age 
of the SNR t, the AM density p am , and the amount of 
collisionless heating at the reverse shock f3 (defined as 
the ratio of electron to ion specific internal energies at 
the shock front). Among these parameters, we have not 
explored variations in the last one (/3), which tend to en- 
hance the Fe emission, and particularly the flux in the Fc 
K blend, due to the higher electron temperatures brought 
on by collisionless heating. In the case of G337.2— 0.7, 
this effect is disfavored by the conspicuous absence of Fe 
K in the X-ray spectrum. For a discussion o n collision- 
less el ectron heating at SNR shock fronts, see Rakowski 
(2005), for its impact on the ejecta em ission, see sections 
2.2 and 2.4 in IBadenes et alJ lJ2005ajl . 

One signifi c ant advantage of the spectra of 
IBadenes et al.l l)2003l) is that they have an under- 
lying physical model for the dynamics of the SNR. 
However, the underlying model is calculated using ap- 
proximations that need to be taken into account for the 
comparison between the data and the synthetic spectra. 
One such approximation is the fact that the models are 
calculated in one dimension. This is clearly a simplifica- 
tion of a very complex situation, because deviations from 
spherical symmetry are expected both on theoretical 
grounds (the contact discontinuity between the forward 
and reverse shocks is subject to Rayleigh- Taylor insta- 
bilities) and due to the large scale structure observed 
in the X-ray images of G337.2— 0.7 (see sections 3 and 
l(j.3|l Thus, the constraints obtained by the comparison 
between the spatially integrated spectrum and the one 
dimensional models of [Badenes et al.1 l)2003fl are focused 
on the emission measure averaged emission of the SNR 
(i.e., they tend to reflect the properties of the brightest 
regions). Local deviations in the properties of the ejecta 
and blast wave emission from the best-fit synthetic 
spectrum are expected, but the average properties of the 
plasma should be similar to those of the most successful 
models. Since we consider the emission from all the 
elements at the same time, it is unlikely that this kind 
of local deviations would conspire to affect our choice 
of models, but this possibility cannot be ruled out 



completely. Another important approximation in the 
models is that the ambient medium is assumed to have 
a uniform density. Again, this is just a simplification of 
a complex scenario. Theoretical studies indicate that 
the circumstellar medium around a Type la SN should 
be strongly modified by the pre-supernova evolution of 
its bin ary progenitor, which w ould sculpt a bubble-like 
cavity l)Badenes fc Brava 1200 1(1 . This should affec t the 
dynamics of the subsequent SNR (Dwarkadas 2000), but 
the observations of prototype Type la SNRs like Tycho 
or SN1006 show no evidence for su ch strongly modified 
dynamics (Badene s fc Bravcl I2001L , see also Badenes 
2004, PhD thesis). On the other hand, t he fact that 
Type la SNe are not detected in the r adio dWeiler et all 
1989), or X-rays (jSchleeel fc Petrel I1993J) seems to 
indicate that the impact of the progenitor evolution on 
the circumstellar medium is limited. Only two Type la 
SNe have shown signs of Ha emiss ion from interaction 
with circu mstellar material 20 02ic l|Hamuv et al.l 2003) 
and 2005gj l)Prieto et al.ll2005l) . In the absence of strong 
evidence to the contrary, we consider the constant 
ambient density hypothesis as a valid approximation, 
and we note that if this were not the case, our detailed 
comparisons might need to be revised. For reference, 
given emission-measure and metallicity considerations, 
the density of the total ambient medium swept-up in 
the ring is about twice that of the diffuse emission (see 
iDwarkadaa I2005L for discussions of how the mass and 
radius of a wind shell might effect the dynamics of the 
SNR). 

The spectra were calculated using a plasma emis- 
sion code maintained by K. Borkowski, whose inputs 
are the electron temperature and ionization structure of 
each ejecta layer as predicted by the SNR model (see 
IBadenes et al.l 120031 for details). An additional Sedov 
component appropriate for the specified age, ambient 
density and explosion energy was used to model the 
blast-wave contribution. The version of the emission 
code used here is missing calculations for Argon as well as 
some of the Li-like lines for other elements that can be 
important for material far from ionization equilibrium. 
(These were included in the NEI code used in the previ- 
ous section.) In particular we believe the disagreement 
in the line centroid for Ca in the upcoming comparisons 
is an artifact of the atomic data used in this NEI emis- 
sion code. However, for the purpose of comparing the 
models with the observed global spectra, over a discrete 
grid in density and age, these emission models were the 
most straightforward to use and readily available to the 
community. Only the normalization and absorbing col- 
umn density are allowed to vary to find the best fit for 
each grid point. 

The x 2 statistic is inappropriate for these comparisons 
since we are not allowing the parameters to vary freely. 
These models have a more limited range of flexibility to 
fit the observations than the general NEI fits used be- 
fore given that the evolution and abundance pattern are 
entirely constrained. However, we can still use the best- 
fit x 2 value at each grid point as a quantitative guide 
for which models are a closer match. It is also useful to 
record the x 2 value over energy bands within the spec- 
trum to confirm our suspicions based on visual inspec- 
tion of why a certain age and density is "fitting" well 
or poorly. These numbers are not of rigorous statistical 
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value, but rather just a rough measure of which models 
are good matches and why. 

Just like the fits presented in Section 5, the goal is 
to match the abundances, ionization state and temper- 
ature of the observed spectrum. However in this case, 
with two connected but separately evolving components, 
matching the data is a matter of finding the best bal- 
ance of these components for any given model. For a 
given density, as time progresses the Sedov component 
becomes more dominant, reducing the apparent equiva- 
lent widths of the metal lines, and the blast-wave compo- 
nent also cools, eventually to below the temperature we 
measured in the last section. Meanwhile the ejecta com- 
ponent is also evolving. The reverse shock is penetrating 
deeper into the layers of nucleosynthetic products, and 
the already-shocked material is continuing to equilibrate 
in temperature and ionize. 

The first test a model must pass is whether there is 
a time when the abundance pattern of the ejecta (mod- 
ified by the solar abundance contribution for the blast 
wave) matches that inferred from the data and the rel- 
evant emission lines are of appropriate strength relative 
to the continuum (from both the Sedov and ejecta com- 
ponents). Three explosion models can be rejected from 
this consideration alone (see Table 0). A pure detona- 
tion explosion (DET) efficiently burns all material to Fe- 
group elements as it propagates supersonically through 
the white dwarf, leaving only a thin skin of intermediate 
elements before it is quenched in the outer layers. Thus 
for DET the iron produced always dominates the Si, S, 
and Ca present in the spectrum, entirely unlike the lack 
of iron (relative to most la models) found in G337.2— 0.7. 
The closest match occurs once the solar abundance ambi- 
ent medium is dominating (Figure (Hp.) which we already 
know cannot produce strong enough lines to reproduce 
the spectrum. 

In the sub-Chandrasekhar model (SCH) a detonation 
originates on the surface of the star and burns the outer 
layers of accreted material. When the ensuing pressure 
waves converge on the center of the WD, a second ig- 
nition happens and a detonation propagates back out. 
Both the center and the outside of the star are incin- 
erated to Fe-group elements, while the middle region is 
less completely burned. For SCH, there simply isn't a 
time at which the abundance pattern is a good match 
for G337.2— 0.7. Prior to ^1000 years the Sedov compo- 
nent is insufficient to fill in the soft continuum. Around 
~1000 years the Si lines are reasonably well matched in 
ionization state and strength relative to the continuum 
(see Figure IHb) but no other portion of the spectrum is 
well-modeled at this time. After ^1000 years the Si and 
S ejecta are swamped by the Sedov component. This 
blast-wave dominated nature occurs earlier for the SCH 
models than it does for the other Type la models because 
of the low ejected mass and explosion energy. 

The pure deflagration models (DEF) are also a poor 
match to the spectrum of G337.2— 0.7. Here the flame 
front is sub-sonic, and the the star expands in reaction 
to the flame propagating out from the center. The flame 
is quenched in the outer regions when the expansion of 
the material is of the same order as the flame velocity. 
This front burns the interior efficiently to Fe-group ele- 
ments but leaves a large portion of unburnt C and O in 
the outer layers. The high column density to G337.2— 0.7 
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Fig. 8. — Type la SNR model comparisons over the original 
grid of densities and ages. The blast-wave (red) and ejecta (blue) 
both contribute to the overall model spectrum (black). The clos- 
est matches for each explosion type are plotted here against just 
the MOS1 data for clarity. From top to bottom, DET pure det- 
onation, pam. = 1.0 X 10 -24 g cm -3 , age = 5000 years; SCH 
sub-Chandrasekhar, pam = 1.0 X 10 -24 g cm -3 , age = 1000 
years. DEFa pure deflagration, p am = 10 X 10 — g cm -3 , age 
= 5000 years. DDTe, delayed detonation, p am = 1.0 X 10 -24 g 
cm -3 , age = 5000 years and PDDe, pulsed delayed detonation, 
Pam = 1.0 X 10~ 24 g cm -3 , age = 5000 years. The bottom two 
panels show the delayed detonation models that are a better match 
to the continuum shape, line fluxes and ionization state than the 
others. 
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masks the emission lines of carbon and oxygen so that 
the preponderance of this unburnt material in the pure 
deflagration models manifests itself as strong continuum 
emission from 1 to 8 keV. Thus even within the ejecta 
component the equivalent widths of Si, S, and Ca are 
fairly low. The Si-rich material lies deep enough in the 
star that by the time Si emission is prominent the Sedov 
contribution is starting to dominate (Figure |5fc) . Fur- 
thermore since it has been shocked relatively recently the 
ionization state of Si is not sufficiently advanced. Soon 
after the DEF "best match" the Sedov contribution dom- 
inates and the Fe-K emission becomes too prominent as 
well. Simulations with a finer grid in ambient density 
and age for the SCH and DEF models did not improve 
upon the fits shown in Figure [S| 

The remaining viable models are the delayed detona- 
tion and pulsed delayed detonation models, DDT and 
PDD, that were preferred initially because of the likeness 
of their overall abundance pattern to that of G337.2— 0.7. 
In delayed detonation models, the flame begins at the 
center as a deflagration, briefly forcing the expansion of 
the star, but at some prescribed flame density it transi- 
tions to a detonation front. Like the DEF models the cen- 
ter is efficiently burned to Fe-group elements, but due to 
the initial expansion the detonation burning in the outer 
layers encounters a wider region where intermediate mass 
elements can be produced before the flame is quenched. 
Thus there is a shell of Si-rich material outside the Fe- 
rich core. The PDD models invoke a pulsation of the star 
to induce the transition from deflagration to detonation. 
They produce a similar pattern of elemental production 
to DDT with a modified density structure. In both cases 
a range of transition densities were explored, with the 
models labeled from a to e in order of the earliest to the 
latest transitions from deflagration to detonation. 

As seen in the previous section, the primary dis- 
agreement between the overall elemental production and 
the inferred abundances in G337.2— 0.7 was the over- 
prediction of Fe. Since the Fe is produced primarily in 
the more interior portion of the star in either delayed det- 
onation model, it is possible that the computed emission 
from Badenes' models might more accurately match the 
G337.2— 0.7 spectrum than the naive comparison to the 
total elemental production. As expected, amongst the 
delayed detonation models, PDDe and DDTe that pro- 
duce the least amount of Fe, are preferred over their more 
energetic cousins when comparing their predicted spectra 
with that of G337.2-0.7. In fact both PDDe and DDTe 
have fairly good matches to the G337.2— 0.7 spectra at 
an ambient density of 1.0 x 10 24 g cm~ 3 and an age of 
5000 years (Table Figure IHJi and e). In comparison to 
the DET, SCH, and DEF models that have "failed", the 
delayed detonation models produce stronger Si, S, and 
Ca lines at the appropriate ionization stage than any of 
the others at close to the right "abundance" ratios. The 
S lines are quite well matched and Si is under-predicted 
but less so than in DEF or DET. At this stage, the reverse 
shock has actually penetrated fully through the ejecta, 
but the Fe remains "hidden" due to a lower electron tem- 
perature and ionization timescale (see below for details) 
and the dominance of the blast- wave component. 

The next question is whether there is an age and den- 
sity where the apparent ionization state of Si, S, and 
Ca are reached before the Sedov component swamps the 
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Fig. 9. — The "best match" models for DDTe and PDDe over 
a finer grid of densities and ages. As in Figure |3 the blast- wave 
(red) and ejecta (blue) both contribute to the overall spectrum 
(black). For DDTe p am = 1.4 X 1CT 24 g cm -3 , age = 6000 year. 
For PDDe p am = 0.8 X 10~ 24 g cm" 3 , age = 7000 year. Only the 
MOS 1 spectrum, which has the highest spectral resolution around 
the Si line-complex, is shown, although all four instruments were 
considered in the analysis. 



ejecta component. To evaluate this requires examining 
a more densely spaced grid in ambient density and age. 
For PDDe, the p am = 1.0 x 10 24 g cm" 3 model, at 5000 
years appeared to be very close to the correct ionization 
timescale, but possibly under-abundant in Si, S and Ca. 
Hence we looked to slightly lower densities where the Se- 
dov solution would take longer to become dominant, but 
towards longer times so that the model would reach the 
same ionization stage (p am = 0.8, 0.9, 1.0 x 10 24 g cm -3 
, age = 5000, 6000, 7000, 8000 years). For DDTe at an 
ambient density of 1.0 x 10 24 g cm~ 3 and age of 5000 
years, the ionization state of Si and S both appear too 
young. There is no evidence for a H-like Si component 
above the Sedov continuum. Hence we looked to slightly 
higher densities and ages around 5000 years for a bet- 
ter match (p am — 1.2, 1.4, 1.6, 1.8 x 10 24 g cm -3 , age = 
4000, 5000, 6000, 7000 years). For all these models, a 
post-facto test of the importance of cooling shows that 
neglecting radiation is a valid approximation for these 
densities and ages. 

Using the \ 2 value relative to all three XMM-Newton 
spectra (MOS1, MOS2 and PN) as a guide we find a best 
match for PDDe at p am = 0.8 x 10 24 g cm" 3 , age = 7000 
years (% 2 = 2287) and for DDTe at p am = 1.4 x 10 24 g 
cm~ 3 , age = 6000 years (x 2 = 2526) as shown in Figure 
El (the implications of these particular parameters will be 
discussed shortly). Given that the only variables in these 
models are the ambient density, the age and the absorb- 
ing column density, these are remarkably good matches 
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Best 



TABLE 5 
"matches" for each Type Ia model 



model 


Pam 


age 


rf w d^ 


D r b 


7V H 


norm 


D c 


x 2 




1CT 24 g cm -3 


years 


pc 


kpc 


10 22 cm^ 2 


(at lOkpc) 


kpc 


(just MOS2) d 


DET 


1.0 


5000 


12.2 


15.5-18.2 


3.96 


0.73 


11.7 


1247.31 


DEFa 


1.0 


5000 


9.8 


12.5-14.6 


4.12 


1.99 


7.1 


786.44 


SCH 


1.0 


1000 


5.7 


7.3-8.5 


3.17 


1.4 


8.1 


795.15 


PDDe 


1.0 


5000 


11.7 


14.9-17.5 


3.22 


0.44 


15.1 


672.26 


DDTe 


1.0 


5000 


11.5 


14.6-17.2 


3.78 


0.77 


11.4 


701.68 


PDDe 


0.8 


7000 


14.0 


17.8-20.9 


3.75 


0.94 


10.3 


549.33 


DDTe 


1.4 


6000 


11.5 


14.6-17.2 


3.80 


0.53 


13.7 


627.66 



a thc forward shock radius implied by Badcncs's model at that age and density 

the distance to the SNR derived by comparing rj w d with the angular radius of the radio remnant (2.7' to 2.3') 
c the distance to the SNR derived from the normalization, D nor7n — (norm)~ ' 10 kpc 

Given only as a guide, cannot be interpreted in the usual sense 



for the character of the G337.2— 0.7 spectrum. While 
the Si abundance or flux in the He-like Si K-shell lines 
(~1.86 keV) are less than they should be, the H-like Si 
Lya(2.0 keV) and Ly/3(2.37 keV) are well matched by the 
combination of the blast-wave and ejecta components. S 
and Ca are primarily contributed from the ejecta com- 
ponent. Both compare well with the data in flux, and S 
appears to agree in ionization state as well. Ne, Mg and 
the soft continuum arc dominated and well-modeled by 
the blast- wave. The hard end of the continuum is slightly 
under-predicted but no strong Fe-K line from the ejecta 
is found that would disagree with the G337.2— 0.7 spec- 
trum. Comparing to the previous figure, the failures in 
the other models may become clearer. Once the first re- 
quirement was met, i.e. the lack of a prominent Fe-K 
line, the most diagnostic aspects of the spectrum were 
the soft-continuum, H-like Si-K, and the S and Ca abun- 
dances. 

The relative contributions from the blast wave and 
ejecta to individual spectral features are also interest- 
ing to note in light of our previous discussions on the 
interpretation of the global and small region spectra. In 
particular the best-matched delayed detonation models 
attribute essentially all of the Ne and Mg line emission 
to the solar-composition blast wave component, exactly 
in line with our explanation for their abundances in the 
global fits. This is also in keeping with the interpretation 
of the bright ring as a circumstellar interaction given that 
the Ne and Mg emission from region " yellow2" preferred 
solar abundances as opposed to the sub-solar abundances 
relative to Si seen in the remnant as a whole. 

Given that the line fluxes, continuum shape and ap- 
parent abundances here are the sum of multiple compo- 
nents, how do the electron temperatures and ionization 
timescales compare with our naive single NEI model? 
Table lists the temperatures and ionization timescales 
for the Sedov blast-wave component and selected ele- 
ments at the specific age and ambient density that best- 
match the spectrum of G337.2— 0.7. The parameters 
listed for the ejecta are the EM weighted averages of the 
electron temperatures and ionization timescales in the 
regions containing those elements. Given the evolution- 
ary state implied by the ambient density and ages found 
here it is not surprising that fits to the global spectrum 
reflect the blast-wave more than the ejecta. The temper- 
atures and timescales we measured under the assump- 



TABLE 6 

Electron Temperatures and Ionization 

Timescales 

for Two of the Type Ia Models 





Sedov a 




Ejecta b 








Si 


Ca 


Fe 


PDDe c 










feT e (keV) 


0.69 


1.68 


1.83 


1.60 


log(n c t) log(cm — 3 s) 


11.58 


11.12 


10.63 


10.26 


DDTe d 










kT c (keV) 


0.65 


2.11 


2.11 


1.82 


log(n £) log(cm — 3 s) 


11.76 


10.53 


10.53 


10.27 



a Mcan shock temperature and ionization age 
Volume emission measure weighted averages for each el- 
ement 
c p<™: 0.8xl0~ 24 gem" 3 age: 7000 years 
d p am - 1.4xl0~ 24 g cm -3 age: 6000 years 



tion of a single kT e and n e t most closely approximate 
the Sedov component of the Type Ia SNR models, but 
may have been pushed towards higher values of kT e and 
n e t by the bright hot ejecta. Naively one might worry 
that the stratified abundance pattern that minimizes the 
Fe-K line flux would also imply that the Ca should be 
under significantly different plasma conditions than Si. 
However, the EM averaged temperatures and ionization 
timescales for Ca and Si are actually quite similar to each 
other in PDDe, and almost identical in DDTe because 
these elements are largely co-spatial in the ejecta (see Ta- 
ble |H| . On the other hand, the EM averaged quantities 
for Fe are quite different from Si, especially given that 
both the lower electron temperature and the shorter ion- 
ization timescale imply a lower level of ionization. These 
differences combined with the relative dominance of the 
blast-wave component allows the Fe to remain hidden. 
None of the ionization timescales in these Type Ia SNR 
models reach the value from the single temperature NEI 
model, however the Sedov ionization timescales are quite 
similar to the values found in the fainter small regions 
outside the bright ring. 

Badenes's Type Ia SNR models naturally provide two 
determinations of the distance to the SNR by predicting 
both the radius of the outer-blast- wave, which can be 
compared to the angular size of the radio SNR, and the 
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EM normalized to a nominal distance of 10 kpc. Compar- 
ing the two distance estimates for a single Type la model 
(Table |SJ we see that they are at least grossly compara- 
ble with each other (not surprising given that the angular 
size used is from the faint material, while the normaliza- 
tion is largely due to the bright material). The distances 
estimated here for the PDDe and DDTe models from the 
normalization and radius are larger than but similar to 
the upper limit on the distance from the H I absorption 
measurements. Also note that the ages and densities of 
the best-matched Type la models are roughly the same 
as those estimated in Section 15.81 for the far end of the 
allowed distance range. However, unlike the estimates 
in Section 15.31 which only considered emission from the 
blast-wave, here the ejecta component is explicitly con- 
sidered and a pre-defined explosion energy has been set. 
In fact these ages and ambient densities were primarily 
chosen to match the relative importance of Sedov and 
ejecta components seen in the spectrum. Given a rea- 
sonable distance estimate of 10 kpc the radius of the 
bright ring would then be 4.6 pc. While this is consider- 
ably larger than a typical planetary nebula, it is actually 
quite similar to the size of some of t he largest plane- 
tary nebula with measured parallaxes (Borkowsk i et al.l 
1990), consistent with the interpretation of the ring as an 
interaction with a structure in the circumstellar medium. 
Simply increasing the ambient density to increase the 
ionization timescale and decrease the radius will not solve 
the remaining discrepancy in distance because an in- 
creased density also results in an increased luminosity 
which in turn impli es a larger distance to the SNR (ex- 
plained in detail bv lBadenes et al]|20 05b). However, the 
history of the remnant is not likely to be as simple as we 
have modeled. Radiative losses to the shock, energy lost 
from the escape of cosmic-rays, the blast-wave running 
into a dense circumstellar ring of material or ambient 
cloud, would all be consistent with a decreased inferred 
shock energy, a higher timescale and a smaller radius. It 
is likely that more than one of these has occurred given 
that roughly three quarters of the explosion energy would 
need to be dissipated to explain the discrepancy between 
the Type la models and the H I distance range. If on 
the other hand general absorption in the region around 
G337.2— 0.7 has obscured a true absorption feature at 
— 116km s _1 then no modification is required and we are 
left with a best estimate of the distance from the X-ray 
spectrum and radio extent alone. 

8. CONCLUSIONS 

The line-rich thermal ASCA X-ray spectrum of SNR 
G337. 2— 0.7 identified it as an "ejecta-dominated" rem- 
nant (Rakowski et al. 2001). In this paper we examine 
the new XMM-Newton, Chandra, and ATCA observa- 
tions of G337.2— 0.7 to learn as much as possible about 
its progenitor system. 

The radio observations provide an image of the full ex- 
tent of the remnant. At both 1.4 and 5 GHz the remnant 
exhibits diffuse emission throughout a shell 2.3' x 2.7' in 
radius that is limb-brightened in some places but not in 
others. In the interior lies a bright elliptical ring of radio 
emission that must indicate a local increase in either den- 
sity or magnetic field (or both). Compared to the X-rays, 
wc sec at least one outer-rim radio filament with X-ray 
emission, and there is bright X-ray emission related to 



the radio ring but brightest in different places. Radio 
H I absorption measurements also provide a constraint 
on the distance of between 2.0 ± 0.5 kpc and 9.3 ± 0.3 
kpc. 

The X-ray spectrum of SNR G337.2-0.7 has been 
studied in detail both globally and locally. Starting 
from a single temperature NEI planar shock model for 
the spectrum of the whole remnant, the previously seen 
super-solar abundances of Si and S are confirmed and 
further constrained (3.50^Qg2 and 5.04lgg2 times their 
solar abundances, respectively). The most striking new 
features of XMM-Newton and Chandra spectra are 
a strong Ca K-shell line, and the lack of any signifi- 
cant K-shell line from Fe. The lack of an Fe-K line 
implies that the soft X-ray spectrum is coming from 
a cooler bremsstrahlung continuum plus a Mg-K line, 
rat her than unreso l ved F e-L lines as previously argued 
by Rakow ski et al.l l)2001|) . The strong Ca line indicates 
the presence of Ca-rich ejecta, yielding a Ca abundance 
of almost 12 times the solar value in the single temper- 
ature, single timescale model. Moreover the global spec- 
trum requires a [Ca/Si] ratio of at least 2.65 times solar 
at the 99% confidence level. In contrast, the [Ar/Si] ratio 
is found to be sub-solar. 

We looked for but could not find any evidence to doubt 
the abundance pattern of Si, S, Ar, and Ca found in the 
single temperature model. Spectra from small regions 
around the remnant showed variations in electron tem- 
perature and ionization timescale but no clear evidence 
for variations in the abundance pattern. The abundances 
in the bright radio ring are if anything lower than the sur- 
rounding areas indicating an origin in the circumstellar 
medium. For the global spectrum two additional simple 
two-temperature models were explored, one with abun- 
dances tied across components and one that had a solar 
composition plasma in addition to the freely varying one. 
These differed from the single temperature model in the 
relative abundances of Ne, Mg or Fe but not of Si, S, Ar, 
and Ca. 

The abundance pattern seen in the global spectrum 
strongly favors a Type la origin. The composition and 
stratification of the ejecta in delayed-detonation Type la 
models can explain the overall abundance pattern if the 
contribution from a solar composition blast-wave into the 
ISM is included and the reverse shock has only recently 
reached the bulk of the Fe-group ejecta. On the other 
hand, typical core-collapse explosions produce Si-group 
elements in the same nuclear burning processes and at 
approximately their solar ratios, inconsistent with the 
abundances of our single temperature global SNR model. 
Core-collapse explosion models that do have a highly 
super-solar [Ca/Si] also have super-solar [Ar/Si] ratios, 
the opposite of which is seen in G337.2— 0.7. The only 
remaining viable explanation for a core-collapse model is 
if a significant amount of hot Ca-rich material is segre- 
gated from the bulk of the Si-rich plasma. Our exami- 
nation of small regions throughout the remnant cannot 
exclude this possibility especially if the two plasmas were 
unresolved but no spatially distinct hot region with a sig- 
nificantly high [Ca/Si] ratio was found. 

As emphasized above, a single temperature model is 
inherently overly simplistic in that it does not reflect the 
variations in composition and conditions within a SNR 
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either between the blast-wave and the reverse-shocked 
ejecta or within the ejecta itself. Given that our ini- 
tial results support a n origin in a Type la SNe and that 
iBadenese t al. (200JJ) have produced predicted spectra 
from Type la SNRs, we have the opportunity to test 
whether our conclusions hold when actual conditions are 
taken into account. Simulated spectra from a variety of 
Type la explosion models were compared with the global 
spectrum of G337.2— 0.7 to find the ambient density and 
age for each model that best matches the spectrum, i.e. 
that avoids a strong Fe-K line, reproduces the strength 
of the He-like Si-K lines relative to the continuum, ex- 
hibits bright S and Ca lines, and reflects the line shapes 
and line-ratios within a given element. In the end we 
saw that indeed Type la models that did not initially 
match the fitted abundances in total metal production 
also did not match the spectrum at any special time in 
their evolution. The preferred delayed detonation and 
pulsed delayed detonation models could in fact repro- 
duce the generic features of the G337.2— 0.7 spectrum. 
In these simulated spectra the reverse shock has over-run 
the Fe-rich material, but there is still no strong Fe-K line 
because Fe is at a lower EM-weighted average tempera- 
ture than the other elements and the high energy X-ray 
band is dominated by blast-wave continuum emission. 

From the above analyses combined with distance con- 
straints from the radio data, we are also able to infer the 
current evolutionary state of G337.2— 0.7. Simple Sedov 
estimates from the single temperature, single timescale 
fits yield ages of betwe en 750 and 3500 year s. For the 
Type la SNR models of Badcncs ct al. (2003), based on 
the relative contributions of the interstellar material and 
the ejecta, G337.2— 0.7 is about 5000 years old. Both the 
age and normalization of the Type la SNR models im- 
ply that G337.2— 0.7 is probably toward the upper limit 
of its allowed distance range. Hence, while G337.2— 0.7 
may not be "ejecta-dominated" in a literal sense, even 
after 5000 years and with more than 50 M of swept 
up material, G337.2— 0.7 still exhibits clear evidence of 
ejecta. 

To reiterate the main purpose of the paper, analysis 
of the X-ray spectrum of G337.2— 0.7 have allowed us to 



identify it as a Type la remnant and, if our identifica- 
tion is correct, show that low-energy (Fe-poor) delayed 
detonations or pulsating delayed detonations are the pre- 
ferred explosion mechanisms. The conspicuous absence 
of Fe K emission in the X-ray spectrum is a strong indica- 
tion that little Fe is found in the outermost layers of SN 
ejecta, where it would be revealed by the high densities 
and electron temperatures. This precludes any Type la 
SN model with plumes of Fe ( 56 Ni) at high velocities, in- 
cluding sub-Chandrasekhar explosions and three dimen- 
sional deflagrations with well-mixed ejecta. We also find 
that the bright ring of radio emission is likely to be ev- 
idence of the impact of the Type la progenitor system 
on the circumstellar medium. In the case of G337.2— 0.7, 
the high [Ca/Si] ratio, combined with a low [Ar/Si] ra- 
tio found in the single ionization, single electron tem- 
perature models, did indeed point to a Type la origin - 
a scenario which then was shown to be viable under a 
more detailed scrutiny. However, the low abundance of 
Fe, which one might have expected to have come from a 
core-collapse explosion, was actually a red-herring, and 
could be explained in the Type la scenarios once actual 
conditions within the ejecta were taken into account. 
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